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Abstract 

We investigate the viability of the quasi-temporal gauge on the lattice. 
This is a complete gauge fixing condition that can be implemented on the 
lattice at a very low computational cost. As a test case, using the Clover 
action, we have evaluated the (gauge invariant) renormalisation constant of the 
non-singlet axial current, using Ward identities extracted from quark states. 
Our result is in reasonable but not complete agreement with previous values 
obtained from Ward identities both on hadronic states and on quark states in 
the Landau gauge. We observe large fluctuations due to lattice Gribov copies. 
The influence of finite volume effects is expected to be non-negligible in the 
case we are considering. 
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1 Introduction 



Gauge fixing, although in principle not necessary for the calculation of gauge in- 
variant quantities on the lattice, is used in practice in several cases. Moreover, it is 
essential to the computation of gauge dependent quantities. 

Typically Landau and Coulomb gauges are used. These are implemented through 
the iterative minimization of a suitable functional. In this way the gauge condition 
can be satisfied to the required precision and in many cases high accuracy is unnec- 
essary. For example, given the gauge condition ||<9fcj4fc|| 2 = 0, where k = 1, . . .,3 
or k = 1, . . . , 4 for Coulomb and Landau gauge respectively, a precision of about 
10 -4 — 10 -5 is adequate when computing fs in the static approximation with 
smeared correlation functions in the Coulomb gauge [1]. On the contrary, high 
precision (at least 10 -10 ) is required for Landau gauge calculations of gluon and 
quark propagators [2], for identifying lattice Gribov copies [3] or for measuring the 
3-gluon vertex function, in order to extract the running QCD coupling [4]. In these 
cases gauge fixing becomes a time consuming part of the computation, comparable 
to the calculation of quark propagators. 

The lattice quasi-temporal gauge (QT gauge), first proposed in ref.[5] and subse- 
quently formulated rigorously in [6], could be a low-cost alternative 1 . It is defined 
by fixing the Coulomb condition at a given arbitrary time t = to and the temporal 
gauge at all points. It is a complete gauge fixing which, compared to the Coulomb 
gauge, is about T times cheaper to implement on the lattice (T is the length of the 
lattice in the time direction). A possible drawback of this gauge condition is that 
it is not invariant under time translations nor is free from the Gribov ambiguity. 
Of course the last problem is not peculiar to this gauge: Gribov copies occur in 
Coulomb and Landau gauges as well. 

The aim of the present paper is to investigate whether the QT gauge can be 
used reliably in actual lattice computations. For this purpose we have calculated in 
this gauge the renormalisation constant of the non-singlet axial current Za , in the 
spirit of [8], by using Ward Identities (WFs) on quark states. Za is suitable for our 
study because it is a gauge invariant quantity that has already been computed in 
many different ways: 

1. from 1-loop lattice perturbation theory [9, 10]; 

2. from WFs on hadron states in a gauge invariant way [11, 12]; 

3. from WFs on quark states in the Landau gauge [12, 13]; 

4. with a non-perturbative method based on amputated quark Green functions 
in the Landau gauge [14]. 

For the third of the aforementioned determinations, a detailed study of the 
influence of Landau gauge Gribov copies on the measurement of Za was performed 
in [13]. As our method for computing Za coincides with the one used in [12, 13], 
except for changing the gauge condition from Landau to QT, we have decided to 
use the same ensemble of configurations that was used in these works, in order to 
allow a detailed comparison of the results. Also, we have analysed the role of the 
lattice Gribov copies in the QT gauge. 

1 Another interesting gauge, proposed in [7], can in principle also be considered. 
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2 The Temporal and Quasi-temporal Gauges 



In this section we discuss the general properties of the QT gauge. First we look 
into the temporal gauge [15], defined by the condition 

A (x) = Mx, (1) 

where Ao{x) is the gauge field. In the naive path integral formulation Gauss's 
law is lost; eq.(l) is an incomplete gauge fixing, as it still allows time independent 
gauge transformations. In other words, a gauge transformation Q,(x) such that the 
gauge transformed field A^(x) satisfies eq.(l) is defined up to an arbitrary gauge 
transformation Q,(x,to) at an initial time to'. 



^/ ^ ~~ iq f drAo(x,r) , ^ 

n(x) = Ve Jt ° 1 ; = tt(x, t )P(x; t ,t), 



(2) 



where V stands for path-ordering and P(x;to,t) = nt=t Uo{x,t) is the Polyakov 
line in the time direction from (x, to) to (x,t). Correspondingly, the tree level naive 
gluon propagator in the continuum 



suffers from a non-physical singularity 1/</q- Its regularization with a principal 
value prescription does not lead, to 0(g 4 ) in perturbation theory, to the correct 
exponentiation of the Wilson loop [16]. One way of solving all these problems at 
once is to integrate on time independent gauge transformations [17]. This operation 
is equivalent to restricting the space of states to those satisfying Gauss's law which 
is then imposed weakly on the physical states. 

As an alternative procedure one can enforce the Coulomb condition at a given 
time t = to [5, 6]: 

d-A(x,t o ) = V(£,* ). (4) 

Eqs.(l) and (4) define the QT gauge. This is a complete gauge fixing, Gauss's law 
is satisfied at t = to and hence at all times and the problematic pole is removed 
from the tree level gluon propagator, which has the form 



2. 



<T(At(x,x )A b m (y,y ))> = j 7^^"« 



a , a e -\<l\\xo-yo\ 



--^y-(ko - yo\-\x - 1 \ - 1 2/o - *o I) 



(5) 



Note that, because of the Coulomb condition at to, the QT gauge breaks transla- 
tional invariance in the time direction. This is why gauge dependent quantities, 
like the propagator of eq.(5), are not invariant under time translations. (Transla- 
tional invariance holds for (5) when the time ordering is either xo < to < yo or 
2/o < to < x .) 

Analytic calculations with the above propagator are hard to perform. However, 
in [6] the correct exponentiation of the Wilson loop to 0(g 4 ) in the QT gauge was 
explicitly demonstrated. 



2 For a different caicuiation of this quantity see aiso [18]. 
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3 The Quasi-temporal Gauge on the Lattice 



We now turn to the implementation of the QT gauge on the lattice. Since at 
the time to one needs to set d^Ak = 0, k = 1, . . .,3, we first review briefly the 
implementation of the Coulomb gauge on the lattice. 

Given a thermalized Monte Carlo configuration U on an L 3 x T lattice with 
periodic boundary conditions, on each timeslice t one can define the functional [19] 



where U^}(n, t) = £l(n, t)Uk(n, t)i2^(n + k, t) represents a gauge transformed SU (3) 
link, i2(n,t) being a gauge transformation matrix. A link configuration U n which 
satisfies the lattice analogue of the Coulomb gauge condition at the time t can be 
obtained by finding a minimum of F(t) with respect to gauge transformations £1. 
More generally, it can be shown that all the extrema of F(t) with respect to £1 
correspond to configurations which satisfy the Coulomb gauge condition d^A^} = 
in discretized form. The existence of more than one minimum is related to the 
presence of Gribov copies. 

It is interesting to note that on the lattice one has in general more classes 
of minima than in the continuum [20, 21]. The issue of distinguishing between 
continuum - like minima and lattice artifacts is certainly an interesting challenge, 
but we shall ignore it in the following. We adopt the pragmatic point of view that 
in a numerical simulation both kinds of minima turn up. 

On the lattice, the minima of F(t) can be found numerically by iteration. The 
typical minimization algorithm sweeps through all lattice sites n at a fixed time t 
and performs local gauge transformations which minimize F(t) with respect to the 
local gauge group. 

Obviously, in order to implement the Coulomb gauge on the entire lattice one 
should repeat the above procedure for all the timeslices, i.e. T times. On the other 
hand, in the QT gauge, we only need to impose the Coulomb gauge on a single 
timeslice to- In this work we have chosen to = 0, i.e. the first timeslice of the 
lattice. In our notation the time index t varies from to (T — 1). 

Once the Coulomb condition holds at t = to, the temporal gauge (1) can be 
trivially imposed by visiting sequentially each timeslice and gauge transforming the 
temporal links Uo(n,t) into the unit group element. On a periodic lattice, this can 
be done for all but one time tj , so that the temporal links Uo(n,tj), rather than 
being unity, end up carrying the value of the Polyakov loops. 

In this work we have chosen tj = T/2 — 1. We have completed the QT gauge by 
fixing first the temporal gauge on the timeslices t = 1, . . . , (T/2 — 1), sweeping the 
lattice in the forward direction, which yields Uo(n, t — l) = l for t = 1, . . . , (T/2 — 1). 
Then the gauge is fixed on the timeslices t = (T — 1), . . . , T/2, sweeping the lattice 
backwards, giving Uo(n, t) = 1 for t = (T — 1), . . . , T/2. Thus, the Polyakov line is 
now between the timeslices (T/2 — 1) and T/2. 

In this way, once the QT gauge is fixed, we end up with a lattice which is 
symmetric in the time direction around t = 0. The discontinuity due to the Polyakov 
line is placed at the greatest possible distance from t = where the signal is killed 
anyway by statistical fluctuations due to the strong damping of Green functions. 

Since the computational cost of fixing the temporal gauge is negligible, it follows 
that the QT gauge is roughly T times faster to implement than the Coulomb gauge. 
This may be a considerable saving in cases where high accuracy is required. 

The QT gauge suffers, like other gauges, from the Gribov ambiguity. Gribov 
copies can be generated when fixing the Coulomb gauge at t = to- When the tempo- 
ral gauge is subsequently implemented, the transformation SI, necessary to perform 
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(6) 



k=l n=l 
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the gauge rotation, will depend for all values oft on possible Gribov copies, so that 
their existence affects the gauge fixing at all timeslices. This is to be contrasted 
to the situation in the Coulomb gauge, where gauge fixing is an independent pro- 
cess on each timeslice and Gribov copies on one timeslice do not depend on what 
happens at other times [3]. 



4 The Measurement of Z 



A 



In order to test the feasibility of lattice non-perturbative calculations in the QT 
gauge, we have calculated the renormalisation constant Za of the axial current 
by using WFs on quark states. As we have pointed out in the introduction, this 
quantity has already been measured with several different methods, using both the 
Wilson and the Clover actions. We adopted the latter action [22, 23] which is free 
from 0(a) discretization errors. We only give an outline of the method and we refer 
the reader to [12] for the details. 

The gauge fixed calculation of Za relies on the following WI for quark Green 
functions 
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2pTr 
pa Tr 



d x dy < u a (y) (u(x)j 5 d(x)) dp(Q) > 



dy< ((j5d(y)) a dp(0) + u(y) a (u(0)j 5 )p) > 



(7) 



In eq.(7) we work explicitly with up (u) and down (d) quark fields with spinor 
indices a and [3. The trace is over colour indices. The value of p is obtained from 
the following ratio of 2-point correlation functions: 



2p: 



d 4 Jdy< A 4 (y,t 9 )vl(d,0) > = 2m 
fdy< V 5 (y,t y )Vl(0,0) > ~ Z A 



(8) 



where 



A^(x) = u{xyf^f 5 d(x) 
V 5 (x) = d(x)j 5 u(x) 



(9) 



are the axial current and the pseudoscalar density and m is the bare quark mass. 
As can be seen from eq.(8), p is gauge invariant. Note that the pa term on the 
r.h.s. of eq.(7) is peculiar to the Clover action and arises from the field redefinitions 
necessary in this formulation [23]. In this work we use the "improved - improved" 
propagators of [10, 24]. 

In order to compute Za we evaluate (8) and the traces in (7) as functions of t y . 
In our case, the (gauge dependent) traces have been evaluated in the QT gauge. 
One can then solve eq.(7) for Za as a function of t y . In the lattice simulations a 
plateau in t y is typically seen, from which the estimate for Za is obtained. 

Actually, barring contact terms, the WI of eq.(7) is expected to hold at all times. 
In practice, as it was already noted in [12], the presence of such terms affects the 
behaviour of the curve near t = 0. 



5 Numerical Results 

We have used data obtained from a Monte Carlo simulation on a 16 3 x 32 lattice 
at [3 = 6.0, with the Clover action of SU{2>) gauge theory, in the quenched ap- 
proximation. An 8-hit Metropolis algorithm was used to generate an ensemble of 
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18 configurations, each separated by 1000 sweeps, after an initial thermalization of 
3000 sweeps. The accuracy of the Coulomb gauge fixing at t = is determined by 
the requirement that F be minimized within a precision of AF/F < 10 -10 , where 
AF is the change in F between two successive gauge fixing sweeps. To increase the 
convergence of the algorithm we used the overrelaxation method [25], setting the 
overrelaxation parameter to u> = 1.72. The whole calculation (and, in particular, 
the gauge fixing) was done in double precision (64 bit). Typically, about 100 gauge 
fixing iterations were needed in order to reach the required precision. In our no- 



1.5 



i i i i i i i i i i i i r 



<> <> <> 



J I I I I I I I I I I I I L 



<> <> 



4 6 8 10 12 14 16 18 20 22 24 26 28 

t 



Figure 1: Za as a function of time in the QT gauge. The value at t = 16 is not 
shown, as it is affected by large fluctuations. 

tation (T = 32 , < t < 31 and tj = 15) the Polyakov line is between f = 15 and 
t = 16. 

In this exploratory study we have used quark propagators at a single value of 
the hopping parameter, k = 0.1425, which corresponds to a pion mass of about 
900 MeV. The value of pa is about 0.05. For light quarks, the dependence of the 
measured renormalisation constants on the quark mass is very mild [24, 26]. 

In order to enhance the signal, we have averaged eq.(7) over four contributions, 
corresponding to the values of the Dirac indices (a, /3)=(1,3), (3,1), (2,4) and (4,2). 
These were found to yield the clearest signal (the same is true in the Landau gauge 
[12]). The statistical errors have been obtained with the jacknife method by deci- 
mating one configuration at a time. 

Our estimate for Za is shown in Fig.l as a function of time t. Except for 
the behaviour near the origin, attributed to the presence of contact terms, in the 
intervals 7 < t < 11 and 22 < t < 26, plateaux appear to settle in. The quality 
of the signal worsens at t ~ T/2 due to the exponential decay of the correlation 
functions. At times t = 15 and 16, where the Polyakov line is located, the signal is 
completely lost. We estimate 

Z A = 1.33 ±0.13 Z A = 1.24 ±0.10 (10) 

from the plateaux at the left and at the right of t ~ T/2 respectively. As our choice 
of imposing the Coulomb condition at to = and the way we fix the temporal gauge 
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both preserve the symmetry of the quark propagators S(x, 0) around t = 0, we can 
average the above results to obtain 

Z A = 1.28 ±0.11. (11) 

Since Za is gauge independent, it is of particular significance to compare this 
result with the Landau gauge one, as they have both been obtained from the same 
WI. In [12] Za was calculated at the same [3 and k values, on the same ensemble 
of 18 configurations. The calculation was performed both with the gauge invariant 
procedure based of WFs on hadron states, and with the gauge dependent procedure 
used in this work, but employing the Landau gauge. The comparison of the values 
obtained for Za as a function of time in the Landau and in the QT gauge is shown 
in Fig. 2 . Although the behaviour of Za in both gauges is qualitatively similar, the 
central values of the QT results are shifted upward with respect to the Landau ones. 
The QT results also display larger statistical fluctuations and a shorter plateau. 
This is reflected in the comparison of the result of eq.(ll) to those quoted in [12]: 
Za = 1-09 ± 0.03 (gauge invariant method) and Za = 1-14 ± 0.08 (Landau gauge). 
We consider the gauge invariant value as the best estimate of Za on this ensemble; 
it is free from the problems associated with gauge fixing and is obtained from a 
much wider and more stable plateau. We note that the QT estimate, although 
compatible with the Landau one, compares rather poorly to the gauge invariant 
one. 
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Figure 2: Comparison of Za calculated in the QT gauge and in the Landau gauge 
[12]. The Landau results were obtained by symmetrising the correlation functions 
with respect to the origin, whereas the QT ones are obtained by averaging the value 
of Za around the origin. 

For completeness, we also quote the results of [13], obtained on a enlarged en- 
semble of 36 configurations at the same lattice volume, [3 and k. From the gauge 
independent WFs they obtain Za = 1-06 ± 0.02, whereas from the WFs on quark 
states in the Landau gauge Za = 1.08±0.05. The latter result has been obtained by 
averaging over 6 Gribov copies for each thermalized configuration. This procedure 
has been introduced in a different framework in ref.[27]. Finally, we also quote from 
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[10] the result Za = 0.97, obtained in boosted lattice perturbation theory in the 
spirit of ref.[28]. 

The discrepancy between the determination of Za in the QT gauge and the 
other estimates reported above could be attributed to the breaking of translational 
invariance in time direction inherent in this gauge. It is conceivable that the loss of 
this symmetry may cause appreciable finite volume effects causing large systematic 
errors in the measurement. 

Another source of error may be the presence of lattice Gribov copies, which, 
as pointed out in [13], also cause an increase of the statistical error. According to 
the discussion in Sect. (3), Gribov copies in the QT gauge affect the gauge trans- 
formations £1 over the whole lattice and may thus have a significant effect on the 
determination of Za ■ 

We have studied the effect of the existence of lattice Gribov copies on Za fol- 
lowing ref.[13]. For each thermalized configuration three lattice Gribov copies have 
been generated by performing random gauge transformations before fixing the QT 
gauge. We have obtained in this way three measurements of Za, using a different 
Gribov copy for each configuration. Besides the result given in eq.(ll), we obtain 
for the three copies: 

Z A = 1.23 ± 0.13 
Z A = 1.30 ± 0.13 

Z A = 1.23 ±0.16. (12) 

Although the results of eqs.(ll) and (12) are compatible, there are big fluctuations 
between results obtained on different copies. Like in [13] our best estimate for Za 
is obtained by averaging over the four values: 

Z A = 1.26 ±0.13. (13) 

This analysis shows that the presence of lattice Gribov copies in the QT gauge is 
responsible of visible fluctuations on the Za value, at the same level of what is 
observed in the Landau case, [13]. Nevertheless the larger error found in the QT 
gauge, compared to the Landau case, cannot be attributed completely to this effect. 

6 Conclusions 

We have investigated the feasibility of lattice computations in the QT gauge. The 
cost of gauge fixing is reduced by approximately a factor of T with respect to 
conventional gauge fixings (Coulomb or Landau). The renormalisation constant of 
the axial current Za has been measured. Our final numerical result agrees only 
roughly with other estimates, within errors. A detailed comparison of our results 
to what is obtained from WFs on hadron states or quark states in the Landau 
gauge shows, in fact, that our numbers are systematically characterised by a higher 
central value and a larger statistical error. We have shown that these features can 
be partly related to the lattice Gribov ambiguity. Also, the breaking of translational 
invariance in the time direction may be responsible for an enhancement of systematic 
errors from finite volume effects. We are currently investigating such effects by 
repeating the calculation at a larger volume. 

We are extremely grateful to G.C.Rossi and M. Testa for an early participation 
to this work and for many illuminating discussions; many useful discussions with 
J.P.Leroy, G.Martinelli and B.Taglienti are also gratefully acknowledged. CP. ac- 
knowledges support from PPARC through an Advanced Fellowship. 
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